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The efTects of cluster diffusion on the submonolayer island density A'' and island-size distribution 
Ns (9) (where Ns is the density of islands of size s at coverage 0) are studied for the case of irreversible 
growth of compact islands on a 2D substrate. In our model, we assume instantaneous coalescence 
of circular islands, while the mobility Ds of an island of size s (where s is the number of particles 
in an island) satisfies Ds ~ s~'^. Results are presented for fi = 1/2 (corresponding to Brownian 
motion), /x = 1 (corresponding to correlated evaporation-condensation), and jx = 3/2 (corresponding 
to cluster diffusion via edge-diffusion), as well as for higher values including n — 2,3, and 6. We also 
compare our results with those obtained in the limit of no cluster mobility (/i — oo). In general, we 
find that the exponents x and x' describing the fiux-dependence of the island and monomer densities 
respectively, vary continuously as a function of p. Similarly, the exponent to describing the flux- 
dependence of the coverage 9m corresponding to the peak island-density also depends continuously 
on /i, although the exponent uj' describing the fiux-dependence of the coverage corresponding to the 
peak monomer density does not. In agreement with theoretical predictions that for point-islands 
with jj, < 1 power-law behavior of the island-size distribution (ISD) is expected, for /x = 1/2 we 
find A'^a ~ s^'^ up to a cross-over island-size Sc- However, the value of the exponent r obtained in 
our simulations (r ~ 4/3) is higher than the point-island prediction r = (3 — m)/2. Similarly, the 
measured value of the exponent corresponding to the dependence of Sc on the average island-size 
S (e.g. Sc ~ 5"") is also significantly higher than the point-island prediction ^ = 2/{fj, + 1). For 
/i < 1, a generalized scaling form for the ISD, Ns{9) — 9/ S^^^'' f{s/ S''), is also proposed, and using 
this form excellent scaling of the entire distribution is found for ^ — 1/2. However, for finite fi > 1 
we find that, due to the competition between two different size-scales, neither the generalized scaling 
form nor the standard scaling form Ns{9) = 9/S^f{s/S) lead to scaling of the entire ISD for finite 
values of the ratio R — Di/F of the monomer diffusion rate to deposition flux. Instead, we find that 
the scaled ISD becomes more sharply peaked with increasing R and coverage. This is in contrast to 
models of epitaxial growth with limited cluster mobility for which good scaling occurs over a wide 
range of coverages. 
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I. INTRODUCTION 

Recently, there has been a lot of interest in under- 
standing the scaling behavior in submonolayer island nu- 
cleation and growthflrl One reason for this is that the 
submonolayer growth regime plays an important role in 
determining the later stages of thin-film growth>i"— Of 
particular interest is the dependence of the total island- 
density N and island-size distribution Ns{6) (where Ns 
is the density of islands of size s at coverage 9 and s is 
the number of monomers in an island) on deposition pa- 
rameters such as the deposition flux F and growth tem- 
perature T. 

One concept that has proven especially useful in stud- 
ies of submonolayer epitaxial growth is that of a critical 
island size^ corresponding to one less than the size of 
the smallest "stable" cluster. For example, if we assume 
that only monomers can diffuse, then in the case of sub- 
monolayer growth of 2D islands on a solid 2D substrate, 
standard nucleation theory^ii predicts that the peak is- 
land density Npk and the monomer density A'^i at fixed 
coverage satisfy. 



where Di h is the monomer hopping rate, i is the critical 
island size, Xi = and x'i = 1 — Xi- We note that in 
the case of irreversible island growth (i — I) this implies 
that xi = 1/3 and Xi =2/3. In addition, it has been 
shown that in the absence of cluster-diffusion and in the 
pre-coalescence regime the island-size distribution (ISD) 
satisfies the scaling form,^'^ 
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Npk ^ {Di,hlF)- 



Ni ^ (Di^h/Fy^^ (1) 



where S is the average island size, and the scaling func- 
tion fi{u) depends on the critical island sizei^ 

However, in some cases (such as in epitaxial growth 
on metal(lll) surfaces) it is also possible for significant 
small cluster diffusion to occuri^SrJ^ In addition, several 
mechanisms for the diffusion of large clusters on solid sur- 
faces have also been proposed J^"— In each case, scaling 
arguments predict that the cluster diffusion coefficient 
Dg decays as a power-law with island-size s (where s is 
the number of particles in a cluster), i.e. Dg ^ s^^. 
In particular, three different limiting cases have been 
consideredi^"— - cluster diffusion due to uncorrelated 
evaporation-condensation {p = 1/2), cluster diffusion 
due to correlated evaporation/condensation {fi — 1), and 
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cluster diffusion due to periphery diffusion (/i = 3/2). We 
note tliat the case ji — 1/2 also corresponds to the Brow- 
nian (Stokes-Einstein) diffusion of compact 2D clusters 
in two-dimensions. 

In order to understand the effects of island diffusion 
on the submonolayer scaling behavior a number of sim- 
ulations have previously been carried out. For example, 
Jensen et aP" have studied the effects of island-diffusion 
with /i = 1 on the percolation coverage for the case of 
irreversible growth without relaxation, corresponding to 
islands with fractal dimension df ~ 1.5. More recently, 
Mulheran and Robbie^ have used a similar model to 
study the dependence of the exponent x on the cluster- 
diffusion exponent ^ for values of /i ranging from to 9. 
They found that for small values of ^ the value of the 
exponent (% ~ 0.45) is significantly larger than the value 
(x — 1/3) expected in the absence of cluster diffusion, 
although it decreases with increasing /i. However, the 
scaling of the ISD was not studied. 

Motivated in part by these simulations, Krapivsky et 
g^j23j24 jj^g^yg carried out an analysis of the scaling behav- 
ior for the case of point-islands, based on the correspond- 
ing mean- field Smoluchowski equations Their analysis 
suggests that due to the large amount of diffusion and 
coalescence in this case, for < 1 the total island den- 
sity saturates (corresponding to "steady-state" behavior) 
while the ISD exhibits power-law behavior of the form, 
Ns ~ s""^, where t = (3 — ^)/2 and the prefactor does 
not depend on coverage.^ This power-law dependence for 
the ISD is predicted to hold up to a critical island-size S'c, 
where Sc ^ and ( = 2/(/i -1-1). In contrast, for ^ > 1 
continuous island evolution is predicted, e.g. the total 
island density does not saturate, and as a result no sim- 
ple power-law behavior is predicted for the ISD. Their 
analysis also indicates that for all values of n, one has 
Xi = x'l = 1/2 with logarithmic corrections. However, it 
should be noted that the point-island approximation is 
typically only valid at extremely low coverages. 

Here we present the results of kinetic Monte Carlo sim- 
ulations of irreversible island growth with cluster diffu- 
sion for the case of compact islands with fractal dimen- 
sion df = 2. Among the primary motivations for this 
work are recent experiments^ on the growth of (com- 
pact) colloidal nanoparticle islands at a liquid- air inter- 
face in which significant cluster diffusion has been ob- 
served. Accordingly, in contrast to much of the previous 
work22i2ii^ our model is an off-lattice model. However, 
our main goal here is not to explain these experiments but 
rather to obtain results which may be used as a reference 
for future work. As already noted, if cluster diffusion 
is due to 2D Brownian motion (as might be expected 
at a fluid-interface) then the value of the exponent /i 
(/i = 1/2) is the same as that expected for uncorrelatcd 



^ The expression t = (3 — ^)/2 has also been derived by Cueille 
and Sira2^ and Camacho™ 



evaporation-condensation. However, we also present re- 
sults for /i = 1 (corresponding to cluster-diffusion due 
to correlated evaporation-condensation), /i = 3/2 (corre- 
sponding to cluster-diffusion due to periphery diffusion) 
as well as for higher values oi fj, {fi — 2, 3, 6 and oo). 

This paper is organized as follows. In Sec. II, we de- 
scribe our model in detail along with the parameters used 
in our simulations, while in Sec. HI we discuss the meth- 
ods we have used to enhance the simulation efficiency. 
In Sec. IV we derive a generalized scaling form for the 
ISD which is appropriate for the case of a power-law ISD 
with T > 1, corresponding to /i < 1. We then present our 
results for the scaling of the island-size distribution and 
island and monomer densities as a function of Dij-i/ F, 
coverage, and /i in Sec. V. Finally, in Sec. VI we discuss 
our results. 



II. MODEL AND SIMULATIONS 

For simplicity we have studied a model of irreversible 
aggregation in which all islands are assumed to be circu- 
lar and rapid island relaxation (perhaps due to periph- 
ery diffusion) is assumed. In particular, in our model 
each island or cluster of size s (where s is the number 
of monomers in a cluster) is represented by a circle with 
area As — irdl/A and diameter dg = dis^/^, where di 
is the monomer diameter. In addition, each cluster of 
size s may diffuse with diffusion rate Dg — Dis^^ where 
Di = Di^h5'^ /A is the monomer diffusion rate, Di h is the 
monomer "hopping rate" , and 5 is the hopping length. 
Similarly, we may write = Ds^h 5^/4 where I?s,h is 
the hopping rate for a cluster of size s. 

In order to take into account deposition, monomers 
are also randomly deposited onto the substrate with rate 
F /d\ per unit time per unit area. Since instantaneous 
coalesce and relaxation is assumed, whenever two clus- 
ters touch or overlap, a new island is formed whose area 
is equal to the sum of the areas of the original clusters, 
and whose center corresponds to the center-of-mass of 
both islands. We note that in some cases a coalescence 
event may lead to overlap of the resulting cluster with 
additional clusters. In this case, coalescence is allowed 
to proceed until there are no more overlaps. In addi- 
tion, if a monomer lands on an existing cluster, then that 
monomer is automatically 'absorbed' by the cluster. 

Thus, at each step of our simulation either a monomer 
is deposited (followed by a check for overlap with any 
clusters) or a cluster is selected for diffusion. If a clus- 
ter is selected for diffusion, then the center of the cluster 
is displaced by a distance 5 in a randomly selected di- 
rection. For computational efficiency, and also because 
it is the smallest length-scale in the problem, in most 
of the results presented here we have assumed 5 — di. 
However, we have also carried out some simulations with 
smaller values {5 = 0.5 di and 5 = 0.25 di) in order to ap- 
proach the continuum limit. As discussed in more detail 
in Sec. VI, our results indicate that the dependence of the 
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island and monomer densities on the hopping distance 5 
is relatively weak. 

We note that besides the exponent /i describing the 
dependence of the cluster diffusion rate on cluster-size, 
the other key parameter in our simulations is the ratio Rh 
of the monomer hopping rate to the monomer deposition 
rate (scaled by the ratio of the hopping length to the 
monomer diameter) e.g.. 



Rh 



Di,h 
F 



(3) 



We note that this definition implies that the dimension- 
less ratio R = Di/Fd\ of the monomer diffusion coeffi- 
cient Di to the deposition flux satisfies, 



R = Rh/^ 



(4) 



Our simulations were carried out assuming a 2D square 
substrate of size L (in units of the monomer diameter 
di) and periodic boundary conditions. In order to avoid 
finite-size effects, the value of L used (L = 4096) was 
relatively large, while our results were averaged over 100 
runs in order to obtain good statistics. In order to de- 
termine the asymptotic dependence of the island density 
on coverage and Rh our simulations were carried out us- 
ing values of R'^ — ARh/ir ranging from 10^ — 10^ up to 
a maximum coverage of 0.3 monolayers (ML). In order 
to study the dependence on /i, simulations were carried 
out for /X = 1/2 (corresponding to Brownian diffusion 
or uncorrelated evaporation-condensation), /u = 1 (cor- 
responding to correlated evaporation-condensation), and 
/i = 3/2 (corresponding to periphery diffusion) as well 
as for higher values = 2, 3, and 6) as well as the case 
/i = (X) corresponding to only monomer diffusion. 

In order to obtain a quantitative understanding of 
the submonolayer growth behavior, we have measured 
a variety of quantities including the monomer density 
-^1 = (i'/4) nxjl? (where n\ is the number of monomers 
in the system) as a function of coverage 0, and the to- 
tal island density N ~ (7i'/4) njl? (where n is the total 
number of islands including monomers in the system). 
In addition, we have also measured the island-size dis- 
tribution Ns{9) where Ng = (7r/4) Us/L"^ corresponds to 
the density of islands of size s. We note that the fac- 
tors of 7r/4 in the definitions above take into account the 
fact that the area of a monomer is (7r/4) d^, and as a 
result the densities defined above all correspond to area 
fractions. Similarly, the coverage 6 — X]s>i corre- 
sponds to the fraction of the total area covered by islands 
(including monomers). 



III. SIMULATION METHODS 

While a simple Monte Carlo approach can be used'^'^ to 
simulate the processes of monomer deposition and cluster 
diffusion such a method can be very inefhcient for large 
values of Rh and small values of /i, since the large range 



of island-sizes and diffusion rates can lead to a low ac- 
ceptance ratio. Accordingly, here we use a kinetic Monte 
Carlo approach. In particular, if we set the deposition 
rate F per unit area df equal to 1, then the total de- 
position rate in the system is while the hopping rate 
for a cluster of size s is given by Rg^h — RhS^^- As 
a result, the total diffusion rate for all clusters is given 
by Rt = ^sRs.h (where is the number of clus- 

ters of size s) while the total rate of deposition onto the 
substrate is i^. The probability Pdep of depositing a 
monomer is then given by. 



dep 



Rt + L^ 



(5) 



while the probability of cluster diffusion is Pdiff — 
1 — Pdep- If cluster diffusion is selected, then a binary 
tree-^^ (whose bottom leaves correspond to the total hop- 
ping rate nsRs,h for each size s) may be used to efficiently 
select with the correct probability which cluster will move 
as well as to efhciently update Rt- However, for large Rh 
and small fi the maximum cluster-size can be larger than 
10^ and as a result the computational overhead associ- 
ated with the binary tree can still be significant. 

Accordingly, we have implemented a variation^ of the 
binary tree approach in which a range of cluster-sizes 
are clustered together into a single 'leaf or bin. In par- 
ticular, to minimize the size of the binary tree, starting 
with island-size s > 3 we have used variable bin-sizes 
such that each bin contains several different cluster sizes 
ranging from a starting value i to a value approximately 
equal to 1.2i. Using this scheme allows us to use a bi- 
nary tree with a maximum of 64 leaves and a rejection 
probability of only 10%. To further decrease the com- 
putational overhead, our binary tree grows dynamically 
from 4 leaves to as many as needed. 

By properly selecting the rates in the binary tree 
and the corresponding acceptance probabilities, one can 
ensure that each diffusion event is selected with the 
proper rate. In particular, if we define the rate of bin 
i as TZb^ = niy.Tlmaxfii, where TZmaxM is the maxi- 
mum cluster-diffusion rate in bin bi (corresponding to 
the smallest cluster-size in the bin) and nfc. is the num- 
ber of islands in the bin, then the sum over all leaves may 
be written. 



(6) 



The probability of attempting a diffusion event is then 
given by. 



P 



diff 



L2 



(7) 



while the probability of selecting bin i is given by Pi — 
TZbjTi^ - Once a bin is selected using the binary tree, a 
specific cluster is then selected randomly from the list of 
all the clusters in that bin. This implies that a cluster 
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of size s will be selected with probability Pg = ns/nti- 
Thus, by assuming an acceptance probability for the se- 
lected cluster-diffusion event given by 

Pace ~ , (8) 

l<-hi 

each diffusion event will be selected with the proper rate. 

Since our simulations are carried out off-lattice, one of 
the most time-consuming processes is the search for over- 
laps every time a cluster is moved. While the simplest 
way to carry out such a search is to check for overlaps 
with all other islands in the system, the search time scales 
as L? ^ and as a result it becomes very time-consuming 
for large systems. Accordingly, we have used a neighbor 
look-up tablc'^^ which contains a list of all other islands 
within a buffer-distance of each island. The search for 
overlaps is then carried out only among the neighbors 
on this list rather than over all the islands in the sys- 
tem. The neighbor list is updated whenever the total 
displacement of any island since the last update is larger 
than half the buffer-distance. 

To speed-up the updates of the neighbor table, we have 
also used a "grid" method^ in which our system is di- 
vided into an Ug by Ug grid of boxes of size Ig — Ljug 
and each cluster can be rapidly assigned to a given box. 
Using this method the search for neighbors only includes 
clusters within an island's box as well as the 8 adjacent 
boxes. As a result, the table update time is reduced to 
Sil? Iv?g instead of l? . To further optimize the speed of 
our simulations, the grid size is varied as the average 
island-size increases. 



IV. GENERALIZED SCALING FORM FOR THE 
ISLAND-SIZE DISTRIBUTION 

As discussed in Sec. I, in both simulations and experi- 
ments on submonolayer epitaxial growth, the island-size 
distribution (ISD) is typically assumed to satisfy the scal- 
ing form given in Eq. [2l However, this scaling form has 
been derived^Ti^ on the assumption that there is only one 
characteristic size-scale S corresponding to the average 
island-size, and that the ISD does not diverge for small 
sjS. However, in our simulations of monomer deposition 
and cluster diffusion and aggregation with < 1, we find 
that the ISD exhibits a well-defined power-law behavior 
for small sj S. In addition, the existence of a shoulder in 
the ISD for large s — implies the existence of a second 
characteristic length-scale which scales as Sc ^ . We 
note that this corresponds to an island size-scale such 
that steady-state behavior breaks down, due to the exis- 
tence of mass-conservation and a finite diffusion length. 

In general one would expect this to lead to a more 
complicated two- variable scaling of the form Ns{9) = 
A g{s/S,s/Sc)- However, if the power-law behavior for 
small s/S is well-defined (and r > 1) then it is possible 
to derive a generalized scaling form involving only one 



variable. In particular, we assume that a scaling form 
for the island-size distribution may be written, 

Ns{0)^A{S,d)f{s/S^) (9) 

In order to determine A{S, 9) note that N — 
Es>iNs = e/S = A{S,e)Y.^^J{s/S)^s. Convert- 
ing to an integral this may be rewritten as 9 / S = 
A{St9) S*^ /i/5c /(") du where u — s/S'^. If we now 
assume that f{u) ~ u^'^ for small u and t > 1, then the 
small-M part of the integral dominates and we obtain, 
A = 9/ S^'^'^^ . This leads to the generalized scaling form, 

= (10) 

We note that a similar scaling form (corresponding to the 
special case C = 1) has previously been derived in Ref.lssI 
for the case of the deposition of spherical droplets with 
dimension D > d on & d-dimcnsional substrate. We also 
note that for C = 1 ^^'^ r = 1 (corresponding to the 
critical value of r) the standard scaling form Eq. [2] is 
obtained. 



V. RESULTS 
A. Stokes-Einstein diffusion (/i = 1/2) 

We first consider the case /Lt = 1/2 corresponding to 
Stokes-Einstein diffusion. Fig. [TJa) shows our results for 
the total cluster density N (including monomers) as well 
as for the monomer density iVi as a function of cover- 
age for three different values of i?^ ranging from 10^ to 
10^. In good agreement with the theoretical prediction 
in Refs. 23 and 24 of "steady-state" behavior for /i < 1, 
wc find that both the monomer density A'^i and total is- 
land density N reach an approximately constant value 
beyond a critical coverage 9m- We note that this cover- 
age decreases with increasing Rh, while the peak island 
and monomer densities also decrease with increasing Rh ■ 

The inset in Fig. \V[h) shows our results for the expo- 
nents X ix — 0.46) and cj (oj ~ 0.38) corresponding to 
the dependence of the peak island density Nm and cover- 
age 9„i on R'f^. In qualitative agreement with the results 
of Mulheran et al2i for fractal islands, the value of x ob- 
tained in our simulations is slightly lower but close to 
1/2. This is also consistent with the prediction^'^'^"' that 
for point-islands x should be equal to 1/2 with logarith- 
mic corrections. Fig.fTJb) shows the corresponding scaled 
island density NR'y^^ as a function of the scaled coverage 
9R'f^^ . As can be seen there is good scaling up to and even 
somewhat beyond the value (9R'f^^ ~ 1) corresponding to 
the peak in the island-density. In contrast, replacing the 
scaled coverage by 9R'f^^ as in Ref. '^T, leads to good scal- 
ing aX 9 — 9jn, but the scaling is significantly worse for 
9 ^ 9„i- Also shown is the scaled monomer density NR'^^ 
(where the peak monomer density scales as Ni^pk ~ R'h~^ 
and the coverage corresponding to the peak monomer 
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density scales as 9i „i ^ R'h ^ ^^'^ j ^ uj' ^ 1/2) as a 
function of the scaled coverage OR'^^^ ■ As for the case of 
the island density, there is good scaling up to and even 
beyond the scaled coverage corresponding to the peak of 
the monomer density. We note that in contrast to the 
exponents x ^-nd the exponent 7 does not appear to 
depend on /i. In particular, we find that for all the values 
of ^ that we have studied, the value of 7 (7 ~ 0.45 — 0.47) 
is close to the value (7 = 1/2) expected in the absence of 
cluster-diffusion. 

We now consJ.der the scaled island-size distribution 
(ISD). In Refs. 23' andlH "steady-state" power-law be- 
havior of the form. 



^1= 1/2 



Ns{e) 



(11) 



where r = (3 — ^)/2 was predicted for /x < 1 for island- 
sizes s << Sc where Sc corresponds to the shoulder in the 
ISD for large s. Similarly, the exponent ^ characterizing 
the scaling of Sc as a function of S (e.g. Sc ^ S'') was 
predicted to satisfy the expression C = 2/(/i+l). We note 
that for /i = 1/2 these expressions imply that r = 5/4 
and C = 4/3. Since N ~ R'f^~^ and S = O/N, one has 
d/S ^ R'f^ ^. Accordingly, Eg . [TT] mav be rewritten as. 



Ns{9) - 9/S 



(12) 



Fig. ^a) shows the ISD scaled using this form. As 
can be seen there is reasonably good scaling for s < Sc, 
although the tail of the distribution does not scale. How- 
ever, the measured value of the exponent r (r ~ 4/3) is 
significantly higher than the predicted value. In addition, 
the measured value of C (C — 3/2) is also significantly 
higher than the predicted value. Fig. ^h) shows the cor- 
responding scaling results obtained using the generalized 
scaling form Eg. [TOl and assuming ^ = 3/2 and r = 4/3. 
We note that this implies that. 



(13) 



As can be seen, in this case both the power-law region for 
small s/S as well as the 'bump' for large s/S scale well 
using this form. We note however, that for the smallest 
clusters (e.g. monomers and dimers) there is poor scaling 
due to deviations from power-law behavior for small s. 



B. Correlated attachment-detachment {fj, — 1) 

We now consider the case fJ, = 1 which corresponds to 
cluster diffusion via correlated attachment-detachment. 
We note that this is the critical value for power-law 
behavior of the ISD (which is expected to occur for 
< ^ < 1) and as a result Krapivsky et al^^ii^ have 
predicted "nested" logarithmic behavior for the island- 
density. Since the simulations are not as computationally 
demanding as for /i — 1/2, in this case we have carried 
out simulations up to 6* = 0.3. Fig.|3Ua) shows our results 
for the total island density N and monomer density A^i 
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FIG. 1: (a) Island and monomer densities and A'^i as a func- 
tion of coverage 9 for = 10^ — 10® and fi = 1/2. (b) Scaled 
densities NR'f^^ and NiR'/^'^ as a function of scaled coverage 
{OR'h^ and 6R'i^^ , respectively). Inset shows dependence of 
peak island density A'^^, and coverage 9m on R'f^. 



as a function of coverage for R'f^ — 10^ — 10^. As can be 
seen, while there is a plateau in the island-density which 
appears to broaden and flatten somewhat with increasing 
R'l^, the plateau is not as flat as for the case fi = 1/2, thus 
indicating deviations from steady-state behavior. As for 
the case fj, = 1/2, a plot of the scaled densities NR'y^ 
{NiR'f^) as a function of scaled coverage OR'^^ {OR'^^) 
shows relatively good scaling up to the coverage corre- 
sponding to the peak island-density, although the value 
of X (x ~ 0.45) is slightly lower than that obtained for 
M=l/2. 

We now consider the island-size distribution. As shown 
in Fig. [3Kb), in this case the ISD does not exhibit a well- 
defined power-law behavior. In particular, on a log-log 
plot the ISD is curved with a slope t ~ 2 for small s and a 
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FIG. 2: (a) Scaled ISD for /i = 1/2 obtained using steady- 
state scaling form Eq. 1121 (b) Scaled ISD obtained using 
generalized scaling form Eq. [10] with r — 4/3 and ( — 3/2. 
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FIG. 3: (a) Island and monomer densities A'^ and A'^i as a func- 
tion of coverage 9 for R'^ = 10^-10^ and = 1. (b) Scaled ISD 
for — 1 using generalized scaling form pO|) with C = 1 a-nd 
T — 2. Results correspond to coverages 6 — 0.025, 0.05 ,0.1, 
= 10^ — 10® and fi = 1. Inset shows corresponding scaling 
results obtained using the standard scaling form 



smaller effective slope (r ~ 1) for large s. Similarly, while 
^ ~ 1 its effective value ranges from 1.03 to 1.1 depending 
on the value of R'/^ and coverage. As a result, neither the 
standard scaling form Eq. [2] nor the generalized scaling 
form Eq. [TU] can be used to scale the entire island-size 
distribution. However, using the generalized scaling form 
pn)) with C ~ 1 and t = 2, we find good scaling for small 
s/S (see Fig. [3fb)), although the ISD does not scale for 
large s/S. On the other hand, if we use the standard 
scaling form ^ (which corresponds to the generalized 
scaling form with C = 1 and t ~ 1, see inset of Fig. [31[b)) 
then the ISD scales for s > Sc but not for smah s. We 
note that this lack of scaling is perhaps not surprising 
since for fj. > I there are two characteristic size-scales S 
and Sc, but no well-defined power-law behavior. 



C. Periphery diffusion (/i — 3/2) 

We now consider the case fi = 3/2 which corresponds 
to cluster diffusion via edge-diffusion. Fig. IH^a) shows 
our results for the total island density N and monomer 
density A^i as a function of coverage for R'f^ = 10^ — 10^. 
As can be seen, while there is a plateau in the island- 
density which appears to broaden with increasing R'f^, 
it is not as flat as for the case fj, — 1, thus indicating 
deviations from steady-state behavior. As for the case 
/i = 1, a plot of the scaled densities NR'^^ {NiR'f/^) as a 
function of scaled coverage OR'i^^ {OR'^//^) shows relatively 
good scaling up to the coverage corresponding to the peak 
island-density. 
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We note that for fj, > 1, Krapivsky et al^=^^ have pre- 
dicted that for point-islands there is a continuous loga- 
rithmic increase in the total island density of the form, 



sin(7r//i) 



\n{eR 



1/2) 



(14) 



However, we find that for fi — 3/2 and higher (not shown) 
scaling plots using this form (e.g. NR^^^ as a function of 
[ln(6'i?i/2)]M/2) provide very poor scaling. In particular, 
since x — 0.45, the scaled peak island-density increases 
with R while the peak position also shifts significantly to 
smaller values. 

We now consider the scaled ISD for /i 3/2. Again in 
this case, it is not possible to scale the entire ISD using 
the average island-size S since there are two characteristic 
size-scales but no well-defined power-law behavior. In 
particular, if we use the generalized scaling form Eq. [10] 
with r = 2 and ^ = 1, then reasonable scaling is only 
obtained for the small-s "tail" corresponding to s/S < 
0.1 (not shown). In addition, as shown in Fig.|4l^b), using 
the standard scaling form Eq. [2] neither the tail nor the 
peak scale. We note that the height and width of the 
"power-law" portion of the ISD decreases with increasing 
and coverage, while the peak near s/S = 1 becomes 
higher and sharper. As a result, the power-law portion 
of the ISD is significantly less important than for smaller 
values of /i. In particular, for R'f^ = 10^ and 6 — 0.1, it 
corresponds to only approximately 10% of the area under 
the curve. 

Fig. [5] shows pictures of the submonolayer morphology 
for R'f^ = 10^ and 9 = 0.1 for ^ = 1/2, 1,3/2, and 2. We 
note that the size-scale M of each picture decreases with 
increasing fi so that approximately the same number of 
islands is visible. As can be seen, in qualitative agree- 
ment with our results, there is a very broad distribution 
of island-sizes for /i = 1/2 while the distribution becomes 
narrower with increasing /i. 
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function of coverage Q for ¥t!y^ = \^ — 10^ and p = 3/2. (b) 
Scaled ISD for = 3/2 using standard scaling form (|10p with 
C = 1 and r = 2^. 



D. Scaling of ISD and densities for n>2 

In order to obtain a better understanding of the de- 
pendence of the island density and ISD on the mobility 
exponent /x, we have also carried out additional simula- 
tions for larger values of /i (/i = 2, 3 and 6) as well as 
in the limit /i = oo in which only monomers can diffuse. 
Fig. ISKa) shows the corresponding results for the scaled 
ISD for /i = 2 using the standard scaling form Eq. [2] for 
different values of the coverage 6 and R'^^. As for the case 
[L — i/2 the ISD does not scale, although the "power- 
law" portion for small s/S is significantly reduced. In- 
stead the peak of the scaled ISD increases with increasing 
coverage and i?^. We also note that for i?^ — 10^ and 
B = 0.1, the peak height is significantly higher than for 
/i = 3/2 while the peak position is closer to s/S ~ 1. 
Similar results for the scaled ISD for /i = 3 are shown 
in Fig. injb), although in this case it tends to sharpen 



more rapidly with increasing R!y^ and coverage. These 
results also suggest that, while the scaled ISD may ap- 
proach a well-defined form (independent of coverage and 
i?Jj) in the asymptotic limit of large the correspond- 
ing scaling function depends on [i. Such a ^-dependence 
is consistent with the dependence of the exponent x and 
x' on \x (see Fig. [5]). 

Fig. [7] shows our results for the scaled ISD for = 6 as 
well as in the limit /i = oo in which only monomers can 
diffuse. Somewhat surprisingly, we find that for /i = 6 
the scaled ISD is significantly broader than for /i = 2 and 
/i = 3, although it is still more sharply-peaked than for 
/i = oo. These results suggest that, at least for (finite) 
fixed , the peak- height depends non-monotonically on 
/i, e.g. it increases from /i = 3/2 to /i = 3 but then 
decreases for higher [i. This is also consistent with our 
results for /I = cx) (see Fig. [71[b)) for which good scaling 
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FIG. 5: Pictures (size M x M) of the submonolayer morphol- 
ogy at coverage 6* = 0.1 and R'^ = 10^ for (a) = 1/2 (M = 
4096) (b) /i = 1 (M = 709) (c) fi = 3/2 (M = 624) (d) = 2 
(M = 485). 



is observed but with a peak height virhich is lower than 
for /X = 6. 

Fig. El^a) shows a summary of our results for the 
monomer density Ni and total island density iV as a 
function of coverage for /i = 1/2,1,3/2,2,6, and oo for 
the case R'/^ — 10^. As can be seen, up to the coverage 
9i^m corresponding to the peak monomer density both 
the island and monomer density are essentially indepen- 
dent of /i. Fig. [HJa) also shows clearly that both the 
island-density and the coverage 9m corresponding to the 
peak island-density increase with increasing n, while the 
monomer density decreases with increasing /i. 

Fig. [Sfb) shows a summary of our results for the de- 
pendence of the exponents x, x'l and 7 on fi. As can be 
seen, the exponent x depends continuously on /i, decreas- 
ing from a value close to 1/2 for small ^ = 1/2) and 
approaching a value close to 1 /3 for large ^. We note that 
these results are similar to previous results obtained for 
fractal islands with df — 1.5 by Mulheran and Robbie^^i 
Similarly, we find that the exponent x' describing the 
dependence of the monomer density at fixed coverage on 
R'f^ also shows a continuous variation with increasing /z, 
starting at a value close to 1/2 for /i = 1/2 and increas- 
ing to a value close to 2/3 for large fi. In contrast, the 
exponent 7 describing the flux-dependence of the peak 
monomer density is close to 1/2 for all /i. 



VI. DISCUSSION 

Motivated in part by recent experiments on colloidal 
nanoparticle island nucleation and growth during droplet 
evaporatioufSS, we have carried out simulations of a sim- 




s/S 

FIG. 6: Scaled ISD for R'^ = 10^ - 10^ and coverage 8 = 
0.025 - 0.1 for (a) ^ = 2 and (b) p = 3. 

plified model of irreversible growth of compact islands 
in the presence of monomer deposition and a power-law 
dependence (Dg ^ s~'^) of the island mobility Dg on 
island-size s. In particular, we have considered the cases 
fi = 1/2 (corresponding to cluster-diffusion via Brownian 
motion), /i = 1 (corresponding to cluster-diffusion via 
correlated evaporation-condensation), and fi = 3/2 (cor- 
responding to cluster-diffusion via periphery diffusion). 
For comparison, we have also carried out simulations for 
higher values of fi including /i = 2, 3 and 6 as well as 
fi = 00. 

In ag reement with the predictions of Ref. [2^ and 
Ref. [2J for point-islands, we find that for small values 
of fi the value of the exponent x characterizing the de- 
pendence of the peak-island density on i?^ is close to but 
slightly lower than 1/2. However, we also find that x de- 
creases continuously with increasing /x, approaching the 
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FIG. 7: Scaled ISD for i?^ = 10^ - 10^ and coverage Q = 
0.025 - 0.1 for (a) ^ = 6 and (b) = oo. 

value 1/3 for large /i. As already noted, these results 
are in good agreement with previous results obtained for 
fractal islandsj^i Similarly, the exponent x' characteriz- 
ing the dependence of the peak monomer density on i?^ 
is also close to 1/2 for small /i, but increases with increas- 
ing fi, approaching the value 2/3 in the limit n ^ oo. In 
contrast, the exponent uj describing the dependence of 
the coverage Om (corresponding to the peak-island den- 
sity) on R'f^ is significantly smaller than 1/2 for small /i 
and also decreases with fi, approaching zero in the limit 
of infinite fi. This is consistent with the fact that when 
only monomers are mobile = oo) the peak island- 
density occurs at a coverage which is independent of R'/^ 
in the asymptotic limit of large R'f^. For comparison, we 
note that while the monomer density Ni{9) depends on 
R'l^ it only depends on /x for coverages beyond the peak 
monomer density (see Fig. [5Ka)). As a result, the ex- 



ponents 7 and w' corresponding to the dependence of 
the peak monomer density (and corresponding coverage 
0i.m) on R'f^ are close to 1/2 for all /i. 

The similarity of our results for x ^-nd to previous 
resultsSi for fractal islands suggests that these exponents 
(along with the exponent x') depend primarily on the 
cluster-mobility exponent fi and substrate-dimension d 
but not on the shape or fractal dimension of the islands. 
We note that such a result is not entirely surprising, 
since for the case in which only monomers can diffuse 
(/Lt = Qo) it has been found that the exponent x depends 
only weakly on the island fractal dimensioni^ In addition, 
we have found that the scaled island and monomer densi- 
ties {NR^ and NiR^) lead to reasonably good scaling as 
a function of OR'"^ , up to and somewhat beyond the peak 
island-density. We note that this scaling form is some- 
what different from that used in Ref. |2l| in which the 
coverage is scaled by so that only the peak scales. 

In addition to the scaling of the island and monomer 
densities, we have also studied the dependence of the 
island-size distribution (ISD) on the cluster-mobility ex- 
ponent /i. In agreement with the predictio n^'^i^^'^^i^'' 
that for point-islands well-defined power-law behavior 
should be observed for /i < 1, for the case /x = 1/2 
we find a broad distribution of island-sizes with a well- 
defined power-law. However, in contrast to the point- 
island prediction that t = (3 — /x)/2 (which implies 
T = 5/4 for /i = 1/2) the value of r obtained in our 
simulations (r ~ 4/3) is somewhat larger. Similarly, the 
value of the exponent {( ~ 3/2) describing the depen- 
dence of the crossover island-size Sc on S for fi — 1/2 
is also significantly larger than the point-island predic- 
tion C — 2/(/i-|- 1) = 4/3. One possible explanation for 
this is that for compact islands the coalescence rate de- 
creases more slowly with increasing island-size than for 
point-islands due to the increase in "aggregation cross- 
section" with increasing island-radius. However, another 
possible explanation is the existence of correlations that 
are not included in the mean-field Smoluchowski equa- 
tions. In particular, we note that in previous work for 
the case of irreversible growth in the absence of cluster 
diffusion (/x = oo), it has been shown^ that there exist 
strong correlations between the size of an island and the 
surrounding capture-zone. 

We note that in contrast to previously studied growth 
models with only limited cluster-diffusion;^^— in which 
there is a single well-defined peak in the ISD correspond- 
ing to the average island-size S, in the presence of sig- 
nificant cluster mobility there are typically two different 
size-scales S and Sc- As a result, in general it is not pos- 
sible to scale the ISD using just the average island-size 
S. However, for the case fi < 1 (corresponding to well- 
defined power-law behavior up to a critical island-size 
Sc) our results confirm that for compact islands the ISD 
exhibits steady-state behavior. As a result, the power- 
law region corresponding to s < Sc can be scaled using 
Eq. [T^l although the large-s "tail" does not scale. Ac- 
cordingly, we have proposed a generalized scaling form 
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FIG. 8: (a) Island density and monomer density A^i as 
function of coverage for 7?^ = 10® and /x = 1/2, 2, 3, 6, and oo. 
(b) Dependence of exponents Xy x'l ^'Hd 7 on the parameter 



FIG. 9: Island density A'^' (not including monomers) and 
monomer density A''i for 5 = d\, S = 0.5di, and continuum 
limit corresponding to 5 = 0. Inset shows dependence of peak 
monomer density A'^i on 6. 



for the ISD, Ns{e) = 6/ S^+^'^ f{s/ S'^) for the case ^ < 1. 
Using this forni, we have obtained excellent scaling for 
the case /x = 1/2. 

In contrast for /i = 1 , there are still two competing size- 
scales S and Sc, but there is no well-defined power-law 
behavior. As a result, no single scaling form can be used 
to scale the entire ISD. However, we find that the value 
of the exponent C (C — 1) is close to that obtained using 
the point-island expression ^ — 2/(fj, + 1). In addition, 
for small s/S the ISD satisfies Ns{0) ^ s~'^''i'f where 
Teff ~ 2. As a result, we find that the small s/S "tail" 
of the ISD can be scaled using the generalized scaling 
form Eq. [10] with r = 2 and C = 1: while the standard 
scaling form Eq.[2]leads to reasonably good scaling of the 
ISD for s> Sc. 

However, for fi > 1 there is no effective power-law be- 
havior and as a result, neither the general scaling form 



Eq. [TU] nor the standard scaling form Eq. [5] lead to good 
scaling of the ISD for finite i?'^. Instead we find that, us- 
ing the standard scaling form Eq.[2]the fraction of islands 
corresponding to small s/S decreases with increasing R'f^ , 
and coverage, while the peak of the scaled ISD increases 
in height and becomes sharper. As a result, the peak po- 
sition shifts to the left with increasing R'/^ and coverage 
and appears to approach 1 for large R'f^. Interestingly, 
this implies, as shown in Figs. Hl^b), [Bl and [71 that for 
/i > 1 the peak of the scaled ISD is even higher than for 
the case of irreversible growth without cluster diffusion 
(/X = 00). However, our results also suggest that, at least 
for fixed coverage and finite (fixed) R'f^, the peak- height 
of the scaled ISD exhibits a non-monotonic dependence 
on fi, since it increases from /i = 3/2to/i = 3 but is 
smaller for fi = 6. 

It is also interesting to compare our results for /i > 1 



11 



with those obtained by Kuipers and Palmer— who stud- 
ied the scaled ISD for the case of fractal islands, assum- 
ing an exponential dependence of the cluster mobility, 
e.g. Ds ~ -Di^* where C < 1- Because of the rapid decay 
of the mobility with increasing cluster-size assumed in 
their model, the resulting scaled island-size distributions 
(using the standard scaling form Eq. [5]) were much closer 
to those obtained for the case of irreversible growth with 
no cluster mobility (e.g. fx = oo) than the results pre- 
sented here. However, for values of ^ which were not too 
small, they also found some evidence of a small island-size 
"tail" , although it was much weaker than found here. 

It is also interesting to consider the applicability of the 
model studied here to recent experiments by Bigioni et 
for the case of colloidal nanoparticle cluster forma- 
tion during drop-drying. We note that in this case, one 
expects that clusters will diffuse on the droplet surface 
via Brownian motion which implies that fi — 1/2. How- 
ever, one also expects that, due to the relatively weak 
Van der Waals attraction between nanoparticles, in this 
case cluster formation may be reversible. Accordingly, 
it would be interesting to carry out additional simula- 
tions for the case of reversible growth corresponding to a 
critical island-size i > 2. 

Finally, we consider the continuum limit of our simu- 
lations. As already mentioned, while our simulations are 
off-lattice, in all of the results presented so far we have 
assumed a hopping length S equal to the monomer diam- 
eter di. We note that this makes our simulations similar 
to previous simulations^"— i^iSiiH with and without clus- 
ter mobility in which a lattice was assumed. However, it 
is also interesting to consider the continuum limit (5—^0. 
In order to do so, we have carried out additional sim- 
ulations with smaller values oi 6 {6 = di/2 and di/A). 
In general, we find that both the monomer density TVi, 



as well as the density A'"' of all clusters not including 
monomers exhibit a weak but linear dependence on the 
hopping length S (see inset of Fig. EKb)) e.g., 

A(5) =A(0)[1 + «(/.) (<5/di)] (15) 

(where X corresponds either to the monomer or island 
density and X(0) corresponds to the continuum limit). 
Accordingly, by performing a linear extrapolation we 
may obtain the corresponding densities in the continuum 
limit. As shown in Fig. [SI for ^ = 1/2 and /i = 3/2 the 
island-density N' depends relatively weakly on the hop- 
ping length, and as a result there is very little difference 
between our results for S = di and the continuum limit. 
In contrast, the monomer density exhibits a somewhat 
stronger dependence on the hopping length 6i. However, 
in general we find a{fi) < 0.1 while the value of a{ii) 
decreases with increasing ^. In particular, in the limit 
fi — oo in which only monomers can diffuse, we find 
a{oo) — 0.01 (0.07) for the island and monomer density 
respectively. These results indicate that in the continuum 
limit the island and monomer densities are only slightly 
lower than in our simulations. Accordingly, we expect 
that in the continuum limit the scaling behavior will not 
be significantly different from the results presented here. 
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